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Table SI: Pairwise estimations for substitution rates (NR data sets) 



non-ARS ARS 



Locus 


Quantile a 


dN 


dS 




dN/dS 


dN>dS d 


dN 


as 


LO 


dN/dS 


dN > dS 


EL A- A 




0.01° 


0.05 


0.36 


0.36 


449(8.38%) 


0.11 


0.07 


1.37 


1.62 


5356(75.99%) 




1 


0.00 


0.01 


0.38 


0.50 


449 


0.03 


0.03 


0.91 


1.25 


877 




2 


0.02 


0.06 


0.33 


0.32 


0 


0.13 


0.07 


1.38 


1.82 


1684 




3 


0.02 


0.05 


0.36 


0.35 


0 


0.13 


0.08 


1.37 


1.66 


1800 




4 


0.02 


0.08 


0.30 


0.30 


0 


0.14 


0.08 


1.56 


1.76 


1160 


HLA-B 




0.01 


0.04 


0.33 


0.30 


121(4.97%) 


0.12 


0.10 


1.27 


1.19 


1637(65.87%) 




1 


0.01 


0.01 


0.47 


0.48 


121 


0.08 


0.07 


1.32 


1.19 


392 




2 


0.01 


0.03 


0.36 


0.36 


0 


0.13 


0.13 


0.97 


0.94 


336 




3 


0.01 


0.04 


0.28 


0.28 


0 


0.14 


0.12 


1.42 


1.20 


419 




4 


0.01 


0.06 


0.24 


0.24 


0 


0.14 


0.09 


1.40 


1.55 


490 


HLA-C 




0.02 


0.05 


0.39 


0.38 


434(7.24%) 


0.07 


0.03 


1.00 


2.53 


4749(79.22%) 




1 


0.00 


0.01 


0.48 


0.52 


434 


0.03 


0.02 


0.89 


1.38 


928 




2 


0.01 


0.04 


0.32 


0.32 


0 


0.07 


0.03 


0.96 


2.82 


1267 




3 


0.02 


0.06 


0.41 


0.42 


0 


0.08 


0.03 


1.11 


2.87 


1282 




4 


0.03 


0.08 


0.37 


0.37 


0 


0.09 


0.03 


1.08 


2.94 


1266 



a, quantiles of divergence (dS non _ARs)- 

b, average pairwise dN/dS. 

c, bold refers to the average pairwise values for each locus. 

d, percentages correspond to the proportion of pairs for which dN > dS in relation to the total number of pairwise comparisons. 
NR, non-recombinant data set 



Table S2: Likelihoods for site models 



Locus 


h 


h 


2{h-h) 


h 


h 


2(k-h) 


HLA-A 


-4361.1(-3537.1) 


-4217.2(-3386.9) 


287.9*(180.5*) 


-4374.7(-3538.6) 


-4218.5(-3386.6) 


312.3*(190.9*) 


HLA-B 


-6969.1(-3537.1) 


-6624.9(-3386.9) 


180.4*(300.4*) 


-6986.5(-3541.6) 


-6632.8(-3386.6) 


190.9*(310*) 


HLA-C 


-3634.1(-3256.9) 


-3562.1(-3257.5) 


144*(-1.2) 


-3635.5(-3257.9) 


-3562.7(-3204.5) 


145.6*(106.8*) 



NOTE.- hand I2, negative log likelihoods for models Ml and M2 (the same applies to models M7 and M8); 
*P < 0.01; d.f.=2 

Numbers in parenthesis refer to the equivalent estimates for NR datasets. 
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Table S3: Site models likelihood convergence for HLA-A (R) 



Model Initial uj Initial k Ml M2 M2wsMl M7 M8 M8usM7 

0.2 2 -4476.54 -4267.42 * -4480.98 -4267.47 * 

0.4 2 -4476.54 -4267.42 * -4483.70 -4266.10 * 

Fequal 0.6 2 -4476.54 -4265.98 * -4479.93 -4265.99 * 

0.8 2 -4476.54 -4266.09 * -4479.90 -4266.11 * 



Flx4 



0.2 2 -4450.86 -4259.73 

0.4 2 -4450.86 -4259.73 

0.6 2 -4450.86 -4259.73 

0.8 2 -4450.86 -4259.73 



-4456.03 -4259. 

-4459.38 -4259. 

-4456.03 -4259. 

-4456.03 -4259. 



F3x4 



02 


2 


-4361.10 


-4217.16 


* 


-4374.67 


-4218.53 


* 




1 


-4361.10 


-4217.16 


* 


-4374.67 


-4218.53 


* 


0.4 


















2 


-4361.10 


-4217.16 


* 


-4374.67 


-4218.53 


* 


0.6 


2 


-4361.10 


-4217.16 


* 


-4374.67 


-4218.53 


* 


0.8 


2 


-4361.10 


-4217.16 


* 


-4374.67 


-4218.53 


* 



0.2 2 -4320.43 -4192.21 * -4336.08 -4194.25 * 

0.4 2 -4320.43 -4192.21 * -4336.08 -4194.25 * 

F61 0.6 2 -4320.43 -4192.21 * -4336.08 -4194.25 * 

0.8 2 -4320.43 -4192.21 * -4336.08 -4194.25 * 

* significant LRT (P < 0.05). 

NOTE.- k, the transition/transversion rate ratio; Fequal, all codons have the same frequency; Flx4, codon frequencies based on 
the frequencies of each nucleotide; F3x4, codon frequencies based on the frequencies of each nucleotide in each codon position; 
F61, each codon frequency is a free parameter. R, data set with recombinants. 
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Table S4: Site models likelihood convergence for HLA-A (NR) 



Model Initial u Initial k Ml M2 M2usMl M7 M8 M28vsM7 

0.2 2 -3701.92 -3491.17 * -3701.97 -3491.30 * 

0.4 2 -3701.92 -3491.17 * -3701.97 -3491.30 * 

Fequal 0.6 2 -3701.92 -3500.97 * -3702.03 -3491.30 * 

0.8 2 -3701.98 -3491.15 * -3701.97 -3491.28 * 



Flx4 



0.2 2 -3668.14 -3464.40 

0.4 2 -3668.14 -3475.11 

0.6 2 -3668.14 -3464.40 

0.8 2 -3668.14 -3475.11 



-3675.15 -3464.41 

-3668.30 -3464.45 

-3668.30 -3464.41 

-3675.15 -3464.41 



F3x4 



02 


2 


-3537.14 


-3386.90 


* 


-3538.56 


-3386.57 


* 


0.4 


1 


-3537.14 


-3386.90 


* 


-3541.58 


-3386.57 


* 


2 


-3537.14 


-3386.90 


* 


-3538.56 


-3386.57 


* 


0.6 


2 


-3537.14 


-3386.90 


* 


-3541.58 


-3386.57 


* 


0.8 


2 


-3537.14 


-3386.90 


* 


-3538.56 


-3386.57 


* 



0.2 2 -3473.08 -3334.39 * -3484.42 -3333.90 * 

0.4 2 -3473.08 -3334.39 * -3474.05 -3333.90 * 

Fcodon 0.6 2 -3473.08 -3334.39 * -3474.05 -3333.90 * 

0.8 2 -3473.08 -3334.39 * -3474.05 -3333.90 * 

* significant LRT (P < 0.05). 

NOTE.- k, the transition/transversion rate ratio; Fequal, all codons have the same frequency; Flx4, codon frequencies based on 
the frequencies of each nucleotide; F3x4, codon frequencies based on the frequencies of each nucleotide in each codon position; 
F61, each codon frequency is a free parameter. NR, data set after removal of recombinants. 



3 



Table S5: Site models likelihood convergence for HLA-B (R) 



Model Initial w Initial k Ml M2 M2usMl M7 M8 M8usM7 

0.2 2 -7228.90 -6760.75 * -7274.97 -6753.2 * 

0.4 2 -7228.98 -6751.69 * -7251.21 -6753.9 * 

Fequal 0.6 2 -7227.86 -6752.46 * -7239.66 -6754.86 * 

0.8 2 -7228.98 -6751.69 * -7240.78 -6755.04 * 



Flx4 



0.2 2 -7222.87 -6767.71 

0.4 2 -7226.59 -6767.08 

0.6 2 -7223.55 -6767.71 

0.8 2 -7227.79 -6767.08 



-7234.06 -6771.06 

-7236.79 -6771.00 

-6767.71 -6771.62 

-7233.72 -6771.99 



F3x4 



02 


2 


-6966.73 


-6623.72 


* 


-6984.45 


-6934.37 


* 




1 


-6970.05 


-6624.89 


* 


-6983.14 


-6632.79 


* 


0.4 


















2 


-6969.08 


-6624.88 


* 


-6986.54 


-6632.79 


* 


0.6 


2 


-6969.88 


-6624.88 


* 


-6981.93 


-6633.68 


* 


0.8 


2 


-6966.54 


-6625.24 




-6983.3 


-6632.79 


* 



0.2 2 -6862.47 -6535.93 * -6876.08 -6547.08 * 

0.4 2 -6860.18 -6539.27 * -6872.78 -6548.52 * 

Fcodon 0.6 2 -6861.53 -6537.12 * -6903.90 -6549.89 * 

0.8 2 -6860.6 -6537.12 * -6905.63 -6543.72 * 

* significant LRT (P < 0.05). 

NOTE.- k, the transition/transversion rate ratio; Fequal, all codons have the same frequency; Flx4, codon frequencies based on 
the frequencies of each nucleotide; F3x4, codon frequencies based on the frequencies of each nucleotide in each codon position; 
F61, each codon frequency is a free parameter. R, data set with recombinants. 
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Table S6: Site models likelihood convergence for HLA-B (NR) 



Model Initial u Initial k Ml M2 M2wsMl M7 M8 M8vsM7 

0.2 2 -3701.98 -3491.17 * -3701.97 -3491.30 * 

0.4 2 -3702.11 -3491.15 * -3702.03 -3491.28 * 

Fequal 0.6 2 -3701.92 -3491.155 * -3702.03 -3491.30 * 

0.8 2 -3701.92 -3491.155 * -3707.27 -3491.28 * 



Flx4 



0.2 
0.4 
0.6 
0.8 



-3668.84 
-3668.14 
-3668.14 
-3668.14 



-3464.4 
-3464.4 
-3464.43 
-3475.11 



-3668.30 -3464.41 

-3668.30 -3464.45 

-3668.30 -3464.41 

-3668.30 -3464.41 



F3x4 



02 


2 


-3537.14 


-3386.90 


* 


-3547.78 


-3386.57 


* 


0.4 


1 


-3537.14 


-3386.90 


* 


-3541.59 


-3386.57 


* 


2 


-3537.14 


-3386.90 


* 


-3541.58 


-3386.57 


* 


0.6 


2 


-3537.14 


-3386.90 


* 


-3541.58 


-3386.57 


* 


0.8 


2 


-3537.14 


-3386.90 


* 


-3541.58 


-3386.57 


* 



0.2 2 -3473.08 -3334.39 * -3475.46 -3333.90 * 

0.4 2 -3473.08 -3334.39 * -3474.05 -3333.90 * 

Fcodon 0.6 2 -3473.08 -3334.39 * -3484.42 -3333.90 * 

0.8 2 -3473.08 -3334.39 * -3475.46 -3333.90 * 

* significant LRT (P < 0.05). 

NOTE.- k, the transition/transversion rate ratio; Fequal, all codons have the same frequency; Flx4, codon frequencies based on 
the frequencies of each nucleotide; F3x4, codon frequencies based on the frequencies of each nucleotide in each codon position; 
F61, each codon frequency is a free parameter. NR, data set after removal of recombinants. 
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Table S7: Site models likelihood convergence for HLA-C (R) 



Model 


Initial to 


Initial k 


Ml 


M2 


WlvsMY 


M7 


M8 


M8usM7 




0.2 


2 


-3697.36 


-3603.67 


* 


-3698.94 


-3603.92 


* 




0.4 


2 


-3697.36 


-3603.67 


* 


-3699.17 


-3605.09 


* 


Fequal 


0.6 


2 


-3697.36 


-3606.20 


* 


-3697.52 


-3767.94 






0.8 


2 


-3697.36 


-3602.37 


* 


-3697.52 


-3602.61 






U.z 


o 

z 


-O/Ul.dO 


-6 i Uz.zo 




-6 i U4.U 


-ODl l.o 


* 




n a 


o 

z 


-O/Ul.OO 


-o / uz.zy 




-o /ui.yo 


-ODll.O 


* 


Flx4 


0.6 


2 


-3702.28 


-3611.76 


* 


-3701.95 


-3703.4 






0.8 


2 


-3701.02 


-3702.29 




-3701.95 


-3611.8 


* 




02 


2 


-3633.59 


-3634.21 




-3636.53 


-3562.82 


* 






1 


-3633.59 


-3634.21 




-3636.07 


-3562.82 


* 




0.4 




















2 


-3634.08 


-3562.14 


* 


-3635.49 


-3562.73 


* 


F3x4 










* 










0.6 


2 


-3634.59 


-3562.77 


-3636.02 


-3562.82 


* 




0.8 


2 


-3634.08 


-3562.14 


* 


-3638.35 


-3636.18 






0.2 


2 


-3572.77 


-3573.43 




-3575.73 


-3506.68 


* 




0.4 


2 


-3572.77 


-3505.98 


* 


-3575.73 


-3506.68 


* 


Fcodon 


0.6 


2 


-3572.92 


-3505.98 


* 


-3575.87 


-3506.68 


* 




0.8 


2 


-3572.77 


-3505.98 


* 


-3575.73 


-3506.68 


* 




* significant LRT (P < 0.05). 













NOTE.- k, the transition/transversion rate ratio; Fequal, all codons have the same frequency; Flx4, codon frequencies based on 
the frequencies of each nucleotide; F3x4, codon frequencies based on the frequencies of each nucleotide in each codon position; 
F61, each codon frequency is a free parameter. R, data set with recombinants. 
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Table S8: Site models likelihood convergence for HLA-C (NR) 



Model Initial u> Initial k Ml M2 M2vsMl M7 M8 M8usM7 

0.2 2 -3333.10 -3378.3 -3335.15 -3265.59 * 

0.4 2 -3333.13 -3265.59 * -3334.89 -3265.59 * 

Fequal 0.6 2 -3335.08 -3267.85 * -3333.17 -3265.59 * 

0.8 2 -3333.10 -3267.85 * -3334.87 -3265.59 * 



0.2 2 -3321.13 -3321.47 -3321.13 -3321.66 

0.4 2 -3321.17 -3254.27 * -3321.25 -3321.43 

Flx4 0.6 2 -3321.25 -3254.25 * -3321.05 -3321.64 

0.8 2 -3321.13 -3321.64 -3321.05 -3321.50 



F3x4 



02 


2 


-3256.94 


-3257.45 




-3258.21 


-3258.36 






1 


-3256.89 


-3204.55 


* 


-3258.16 


-3258.74 




0.4 


















2 


-3256.89 


-3257.46 




-3257.91 


-3204.64 


* 


0.6 


2 


-3256.89 


-3204.49 


* 


-3257.87 


-3258.66 




0.8 


2 


-3256.94 


-3257.26 




-3257.91 


-3204.50 


* 



0.2 2 -3204.57 -3164.9 * -3205.78 -3208.89 

0.4 2 -3204.57 -3154.39 * -3205.78 -3154.55 * 

Fcodon 0.6 2 3204.57 -3205.05 -3205.78 -3206.28 

0.8 2 3204.57 -3154.39 * -3205.78 -3154.55 * 

* significant LRT (P < 0.05). 

NOTE.- k, the transition/transversion rate ratio; Fequal, all codons have the same frequency; Flx4, codon frequencies based on 
the frequencies of each nucleotide; F3x4, codon frequencies based on the frequencies of each nucleotide in each codon position; 
F61, each codon frequency is a free parameter. NR, data set after removal of recombinants. 



Table S9: Classifications of ARS codons 



Reference 


Total number of codons 


Bjorkman et al. (1987) (BJOR) 


57 


Chelvanayagam (1996) (CHEV) 


48 


Yang and Swanson (2002) (YANG) 


26 


BIT 


33 



BJOR and CHEV are ARS classifications. 

YANG is a list of codons with significant ui > 1 in HLA genes. 

BIT is the set of codons with P > 0.95 (M8; see Methods for details). 
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Table S10: Parameter estimates under site models for EL A- A 



Data set Model Parameter estimates Sites with to > 1 

R MO u = 0.96 None 

Ml p Q = 0.80, p = 0.2, lo 0 = 0.05, Wi = 1 Not allowed 

M2 po = 0.68, pi = 0.25, p 2 = 0.07 9F,62Q,63E,66N,70H,76A,77N,81L,82R,97I, 

Lo 0 = 0.08, wi = 1, lo = 9.12 99Y,114R 

M7 p = 0.08, 9 = 0.20 Not allowed 

M8 po = 0.93, p = 0.09, q = 0.16 9F,62Q,63E,66N,70H,76A,77N,81L,82R,95I, 

{pi = 0.07), cj s = 9.44 97I,99Y,114R,116D,151H,152A,156R,163R 

NR M0 to = 0.93 None 

Ml po = 0.76, p = 0.24, w 0 = 0.06, lj 1 = I Not allowed 

M2 po = 0.63, pi = 0.30, p 2 = 0.07 9F,62Q,63E,66N,70H,76A,77N,81L,97I, 

Lo 0 = 0 09, wi = 1, to = 8.54 99Y,114R,151H,152A,156R,163R,167G 

M7 p = 0.02, q = 0.03 Not allowed 

M8 po = 0.93, p = 0.02, q = 0.02 9F,62Q,63E,66N,70H,76A,77N,81L,97I 

(pi = 0.07), w s = 9.11 99Y,114R,151H,152A,156R,163R, 167G 

NOTE.- F3x4, initial k = 2; positively selected sites by the BEB method were identified based on the p > 95% criterion. 
Sites with p > 99% are in bold. po,Pi,P2, <*>o, w li w 5 proportions of sites within each category of ui values (purifying selection, 
neutral evolution, positive selection); peg, shape parameters for the beta distribution of M7; value estimated from data under 
M8. Letters attached to the site positions refer to the aminoacid present in that position in the reference sequence for the locus. 
R, data set with recombinants; NR, data set after removal of recombinants. 
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Table Sll: Parameter estimates under site models for HLA-B 



Data set Model Parameter estimates Sites with to > 1 

R MO uj = 0.77 None 

Ml po = 0.83, p = 0.17, cj 0 = 0.02, lj 1 = 1 Not allowed 

M2 po = 0.796, pi = 0.13, p 2 = 0.07 24S,45E,63N,67Y,70Q,71A,74D,77S, 

Lo 0 = 0.03, wi = 1, w = 7.66 80N,81L,82R,94T,95L,97S,113H, 

114D,116Y,152E,156R,163E,305A 
M7 p = 0.01, g = 0.02 Not allowed 

M8 p 0 = 0.93, p = 0.01, q = 0.05 24S,45E,63N,67Y,69A,70Q,71A,74D, 

(pi = 0.07), w s = 8.098 77S,80N,81L,82R,94T,95L,97S,113H, 

114D,116Y,152E,156R,163E,305A 
NR M0 cj = 0.86 None 

Ml po = 0.83, p = 0.17, cj 0 = 0.02, wi = 1 Not allowed 

M2 po = 0.77, pi = 0.14, p 2 = 0.09 9Y,11S,12V,24S,45E,63N,67Y,70Q,71A, 

Lu 0 = 0.02, co! = 1, w = 9.64 74D,77S,80N,81L,82R,83G,95L,97S,99Y, 

113H,114D,116Y,152E,156R,163E,171Y,305A 
M7 p = 0.01,g = 0.09 Not allowed 

M8 po = 0.91 p = 0.02 q = 0.14 9Y,11S,12V,24S,45E,63N,67Y,69A,70Q,71A,74D, 

(pi = 0.09) oj s = 9.48 77S,80N,81L,82R83G,95L,97S,99Y,113H, 

114D,116Y,152E,156R,163E,167W,171Y,305A 

NOTE.- F3x4; initial k = 2; positively selected sites by the BEB method were identified based on the p > 95% criterion. 
Sites with p > 99% are in bold. po,Pi,P2, <*>o, w ii w 5 proportions of sites within each category of u> values (purifying selection, 
neutral evolution, positive selection); peg, shape parameters for the beta distribution of M7; value estimated from data under 
M8. Letters attached to the site positions refer to the aminoacid present in that position in the reference sequence for the locus. 
R, data set with recombinants; NR, data set after removal of recombinants. 
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Table S12: Parameter estimates under site models for HLA-C 



Data set Model Parameter estimates Sites with lo > 1 

R MO u = 0.73 None. 

Ml po = 0.74, p = 0.26, cj 0 = 0.04, uj 1 = 1 Not allowed 

M2 po = 0.65, pi = 0.32, p 2 = 0.03 9F,24S,73T,77S,80N,99C,116Y,156R,163T 

uj q = 0.03, coi = 1, w = 9.8 
M7 p = 0.01, 9 = 0.02 Not allowed 

M8 po = 0.97, p = 0.005, q = 0.01 9F,24S,73T,77S,80N,99C,116Y,156R,163T 

{pi = 0.03), uj s = 10.47 
NR M0 cj = 0.74 None 

Ml po = 0.69, p = 0.31, oj q = 0.04, wi = 1 Not allowed 

M2 p Q = 0.69, Pl = 0.31, p 2 = 0.00 9F, 24S, 73T, 77S, 80N, 116Y, 156R, 163T 

lu q = 0.04, wi = 1, u> = 999 a 
M7 p = 0.02, g = 0.03 Not allowed 

M8 po = 0.95, p = 0.05, q = 0.006 9F,24S,73T,77S,80N,99C,116Y,156R,163T 

(pi = 0.05),cj s = 8.79 

NOTE.- F3x4, initial k = 2; positively selected sites by the BEB method were identified based on the p > 95% criterion. 
Sites with p > 99% are in bold. po,Pi,P2, <*>o, w ii w 5 proportions of sites within each category of ui values (purifying selection, 
neutral evolution, positive selection); peg, shape parameters for the beta distribution of M7; value estimated from data under 
M8. Letters attached to the site positions refer to the aminoacid present in that position in the reference sequence for the locus. 
R, data set with recombinants; NR, data set after removal of recombinants. 



Table S13: Correlations between pairwise dN/dS ARS and divergence (dS non _ A Rs) 



Locus 


Spearman 


Kendall 


Pearson 


HLA-A 
HLA-B 
HLA-C 


0.12*(0.22*) 
0.20*(0.19*) 
0.26*(0.24*) 


0.08*(0.14*) 
0.14*(0.13*) 
0.18*(0.16*) 


0.17*(0.25*) 
0.20*(0.12*) 
0.20*(0.19*) 



* P < 0.001 obtained by Mantel test. 

NAs handled by casewise deletion. 

In parentheses: values for NR data sets. 
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Table S14: Branch models likelihood convergence for HLA-A (ARS; NR) 



Frequency Model 


Initial uj 


Model 0 (lnl; uj) 


Model 2a (lnl; uj w ; 




Model 2b (lnl; to Up ; 


Mint) 




0.2 


-946.16; to 


= 4.59 


-946.15; 


to w 


= 4.34; 


UJ b 


= 4.88 


-948.55; uj t = 


6.91; uj, 


= 3.77 




0.4 


-946.16; to 


= 4.59 


-946.15; 




= 4.34; 


UJ b 


= 4.88 


-945.92; cu t = 


3.43; uj, 


= 5.95 


Fequal 


0.6 


-946.16; to 


= 4.59 


-946.15; 


U! w 


= 4.34; 


UJ b 


= 4.88 


-945.92; uj t = 


3.43; uj, 


= 5.95 




0.8 


-946.16; to 


= 4.59 


-946.15; 


UJ W 


= 4.34; 


UJ h 


= 4.88 


-945.92; uj t = 


3.43; uj, 


= 5.95 




0.2 


-943.60;w 


= 4.23 


-943.59; 


to w 


= 4.01; 


UJ h 


= 4.47 


-945.30; to t = 


6.25; tOi 


= 3.50 




0.4 


-943.60;u; 


= 4.23 


-943.59; 


to w 


= 4.01; 


tOb 


= 4.47 


-947.898; co t = 


- 6.75; cji 


= 3.82 


Flx4 


0.6 


-943.60; to 


= 4.23 


-943.59; 


UJ W 


= 4.01; 


UJ h 


= 4.47 


-947.898; co t = 


- 6.75; tUi 


= 3.82 




0.8 


-943.60;u; 


= 4.23 


-943.59; 


U) w 


= 4.01; 


UJ b 


= 4.47 


-947.898; uj t = 


6.75; uji 


= 3.82 




0.2 


-935.21; uj 


= 1.84 


-935.18; 


U) w 


= 1.68; 


UJ b 


= 2.03 


-934.96; uj t = 


1.39; uji 


= 2.35 




0.4 


-935.21; to 


= 1.84 


-935.18; 


to w 


= 1.68; 


UJ b 


= 2.03 


-934.96;w t = 


1.39; uj, 


=2.35 


F3x4 


0.6 


-935.21;w 


= 1.84 


-935.18; 


to w 


= 1.68; 


UJ b 


= 2.03 


-934.96; tu t = 


1.39; uj, 


= 2.35 




0.8 


-935.21; to 


= 1.84 


-935.18; 


UJ W 


= 1.68; 


UJ b 


= 2.03 


-934.96; to t = 


1.39; uj, 


= 2.35 




0.2 


-910.02; uj 


= 1.82 


-909.91; 


UJ W 


= 1.59; 


UJ b 


= 2.17 


-910.01; uj t = 


1.74; uj, 


= 1.89 




0.4 


-910.02; to 


= 1.82 


-909.91; 




= 1.59; 


UJ b 


= 2.17 


-910.01; uj t = 


1.74; uj, 


= 1.89 


Fcodon 


0.6 


-910.02; uj 


= 1.82 


-909.91; 


to w 


= 1.59; 


UJ b 


= 2.17 


-910.01; uj t = 


1.74; uj, 


= 1.89 




0.8 


-910.02; uj 


= 1.82 


-909.91; 


UJ W 


= 1.59; 


UJ b 


= 2.17 


-910.01 ; uj t = 


1.74; uj, 


= 1.89 



Model 2a: within/between; Model 2b: terminal/internal; lnl, log likelihood. NR, data set after removal of recombinants. 
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Table S15: Branch models likelihood convergence for HLA-B (ARS; NR) 



Model Initial uj Model 0 (lnl; uj) Model 2a (lnl; uj w ; uj b ) Model 2b (lnl; uj t i P ', Mint) 

0.2 -1159.28; uj = 2.79 -1158.96m,, = 2.08; oj b = 3.26 -1158.77; u t = 1.9; uj, = 3.38 

0.4 -1159.28; uj = 2.79 -1158.96; uj w = 2.08; u b = 3.26 -1158.77; u t = 1.9; = 3.38 

Fequal 0.6 -1159.28; oj = 2.79 -1158.96; co w = 2.08; w b = 3.26 -1158.77; wt = 1.9; oji = 3.38 

0.8 -1159.28; uj = 2.79 -1158.96; w w = 2.08; w 6 = 3.26 -1158.77; u t = 1.9; = 3.38 



0.2 -1159.21; uj = 2.58 -1158.87; uj w = 1.95; uj b = 2.99 -1158.75; u t = 1.79; w 4 = 3.11 

0.4 -1159.21; uj = 2.58 -1158. 93;^ = 1.95; w 6 = 2.99 -1158.87; to t = 1.95; uj, = 2.99 

Flx4 0.6 -1159.16; uj = 2.58 -1158.93; uj w = 1.95; cj 6 = 2.99 -1158.87; uj t = 1.95; uj, = 2.99 

0.8 -1159.16; uj = 2.58 -1158.93;^ = 1.95; uj b = 2.99 -1158.87; uj t = 1.95; w 4 = 2.99 



0.2 -1112.27; uj = 0.99 -1111.91; w w = 0.73; w b = 1.16 -1111.78; cj t = 0.69; uj, = 1.2 

0.4 -1112.27; cj = 0.99 -1111.91; w w = 0.73; w 6 = 1.16 -1111.78; u t = 0.69; oj, = 1.2 

F3x4 0.6 -1112.27; cj = 0.99 -1111.91; w w = 0.73; w 6 = 1.16 -1111.78; uj t = 0.69; ^ = 1.2 

0.8 -1112.27; cj = 0.99 -1111.91; <*;„, = 0.73; oj b = 1.16 -1111.78; w t = 0.69; ^ = 1.2 



0.2 -1088.75; uj = 1.88 -1088.27; uj w = 1.3; w b = 2.29 -1089.39; w t = 1.20; uj, = 2.36 

0.4 -1088.75; uj = 1.88 -1088.27; uj w = 1.3; w b = 2.29 -1089.8; uj t = 1.27; uj, = 2.31 

Fcodon 0.6 -1090.33; uj = 1.88 -1088.27; uj w = 1.3; w b = 2.29 -1088.76; w 4 = 1.20; uj, = 2.35 

0.8 -1088.75; uj = 1.88 -1089.84; uj w = 1.30; cj b = 2.29 -1088.11; uj t = 1.22; uj, = 2.34 



Model 2a: within/between; Model 2b: terminal/internal; lnl, log likelihood. NR, data set after removal of recombinants. 
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Table S16: Branch models likelihood convergence for HLA-C (ARS; NR) 



Model Initials Model 0 (lnl; co) Model 2a (lnl; uj w ; u; b ) Model 2b (lnl; co t ; iUi) 

0.2 -946. 76;cj = 4.47 -946.89; uj w = 2.77; u = 9.6 -944.33; u t = 2.12; Wi = 11.57 

0.4 -946. 76;cj = 4.47 -945.49;^ = 2.77; w b = 9.6 -944.33; u t = 2.12; Wi = 11.57 

Fequal 0.6 -946. 76;cj = 4.47 -946.89; w„ = 2.77; w b = 9.6 -944.33; w t = 2.12; Ui = 11.57 

0.8 -946. 76;cj = 4.47 -946.89; w w = 2.77; u; b = 9.6 -944.33; w t = 2.12; w 4 = 11.57 



0.2 -960.73;cj = 4.04 -959.47; u w = 2.51; uj b = 8.66 -958.44; u t = 1.97;w 4 = 10.29 

0.4 -960.73;w = 4.04 -959.47; u w = 2.51; uj b = 8.66 -959.47; u t = 2.51; = 8.66 

Flx4 0.6 -960.73;cj = 4.04 -959.47; u w = 2.51; u b = 8.66 -959.47; u t = 2.51; lo, = 8.66 

0.8 -960.73;cj = 4.04 -959.47; u w = 2.51; uj b = 8.66 -959.47; u t = 2.51; = 8.66 



0.2 -951. 50;w = 1.89 -950.199;^ = 1.19; co b = 4.14 -947.82;w t = 0.95 ; cj, = 4.91 

0.4 -951. 50;u; = 1.89 -950.199 u w = 1.19; u b = AAA -947.82; uj t = 0.95; a;* = 4.91 

F3x4 0.6 -951. 50;w = 1.89 -950.199 u w = 1.19; w 6 = 4.14 -947.82; w t = 0.95; = 4.91 

0.8 -951. 50;u; = 1.89 -950.199; u w = 1.19; w 6 = 4.14 -947.82; w t = 0.95; = 4.91 



0.2 -900.96;cj = 2.99 -899.56; lu w = 1.8; u b = 6.65 -898.68; w t = 1.46; = 7.67 

0.4 -900.96;cj = 2.99 -899.56 ; cu w = 1.8; uo b = 6.65 -898.68; u t = 1.46; u t = 7.67 

Fcodon 0.6 -900.96;w = 2.99 -899.56; lu w = 1.8; uj b = 6.65 -898.68; u t = 1.46; u t = 7.67 

0.8 -900.96;cj = 2.99 -899.56; u> w = 1.8; iv b = 6.65 -898.68; u t = 1.46; w 4 = 7.67 



Model 2a: within/between; Model 2b: terminal/internal; lnl, log likelihood. NR, data set after removal of recombinants. 
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